Electromagnetic Spin-Orbit Interactions via Scattering 
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The longitudinal components of orthogonal-circularly polarized fields carry a phase singularity 
that changes sign depending on the polarization handedness. The addition of orbital angular mo- 
mentum adds to or cancels this singularity and results in polarization-dependent scattering through 
round and square apertures, which we demonstrate analytically, numerically, and experimentally 
By preparing the incident polarization and arranging the configuration of sub- wavelength apertures, 
we produce shadow-side scattered fields with arbitrary phase vorticity 



Within the last decade, there have been wide obser- 
vations of optical phase singularities in the evanescent 
fields produced by propagating [1 and scattered [2 light. 
These phase singularities indicate where the electric field 
is strictly zero and carry information about the elec- 
tric field Poynting vector and angular momentum [3 . 
A three-dimensional electric field produces three differ- 
ent types of polarization phase singularities [4], the evo- 
lution of which is studied in a rich array of literature 
[5]. Our understanding of phase singularities allows us 
to probe and characterize materials, surfaces, and light 
propagation dynamics, as well as control and manipulate 
microparticles [6 . 

Recently, there have been reports that the locations of 
near-field phase singularities produced by chiral "gam- 
madion" [7] and spiral grating structures [8 depend on 
incident polarization handedness. These phase singular- 
ities are connected to the extraordinary transmission of 
light through sub- wavelength slits [9 , where whirlpool- 
like power flows and singularities in the Poynting vector 
are shown to exist ^01[11J. Moreover, azimuthally and 
radially-polarized vortices, beams with different polar- 
ization singularities, are transmitted through apertures 
with different efficiencies [T2^. Despite numerous mea- 
surements and observations of phase singularities in the 
near field, the polarization-dependent transmission that 
occurs at sub-wavelength structures is not fully under- 
stood, so that light-metal interactions are neither fully 
optimized nor controlled. 

Here, we show that polarization-dependent phases 
at subwavelength-structured materials are concisely de- 
scribed by electromagnetic spin-orbit interactions, where 
"spin" refers to the circular polarization of electric fields, 
and "orbital angular momentum" identifies vortex phase- 
fronts that are the signature of phase singularities. Elec- 
tromagnetic spin angular momentum and orbital angular 
momentum are coupled specifically when the longitudi- 
nal electric field component plays a role in the dynamics 
and previously cited examples of optical spin-orbit inter- 
actions investigate oblique reflection or refraction flU [15] 



or high-numerical aperture focusing [16j. 

We explain, for the first time, that polarization- 
dependent singularities in longitudinal field components 
dictate which modes and to what extent light is trans- 
mitted through thin-film apertures via the process of 
electromagnetic scattering. Our results suggest that 
waveguiding due to finite-thickness materials [12 is not 
the only explanation for understanding polarization- 
dependent transmission through round subwavelength 
apertures. Numerical simulations verify our analyti- 
cal predictions and we reconstruct experimental THz 
field measurements with subwavelength-resolution that 
demonstrate polarization-dependent phase singularities. 
Our new insight of coherent light-metal interactions en- 
ables us to produce shadow-side longitudinal fields with 
arbitrary phase vorticity by controlling the polarization 
and configuration of sub- wavelength apertures. 

The topological charge or phase winding associated 
with the longitudinal or 2:-component of a spin-polarized 
electric field is = mi -\- m^, where mi is the topo- 
logical charge associated with orbital angular momen- 
tum in the transverse electric field components, and 
the topological charge associated with the photon spin 
number, = ±1, depends on the orthogonal-spin po- 
larization a±. The longitudinal component topological 
charge is illustrated by representing a continuous-wave 
circularly-polarized field with transverse mode amplitude 
A(/),0,z)e^^^^ as 
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where the circular polarization unit vector is represented 
in cartesian and cylindrical coordinates as (T± = (i ± 
ij)/\/2 = {p±i4>)e^'^^ /V^. The longitudinal component 
of the electric field is calculated by Maxwell's equation 
V-E = (V^+5,)-E = 0, 
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where = 6>, z^). Equation [s] indicates that an elec- 
tric field distinctly gains or loses unit topological charge 
rris = ±1 when represented in the circular polarization 
basis. The vortex phase coefficient exp[±i0] is inter- 
preted as a path-dependent geometric phase [171 [IH] ? 
it describes the phase that is accumulated around an az- 
imuthal path in cylindrical coordinates, or a star singu- 
larity at a C-point, where the lines depicting pure circular 
polarization curve due to the finite transverse spatial ex- 
tent of the electric field [4 . Spin-orbit interactions arise 
because the spin and orbital angular momentum contri- 
butions, generally associated with the first and second 
round-bracketed terms of Eq. [Sj respectively, add con- 
structively to the longitudinal electric field magnitude. 
Since the orbital angular momentum contribution scales 
inversely with radius, spin-orbit interactions are signifi- 
cant when the distances between phase singularities and 
scattering edges are less than the incident field wave- 
length. 

We use an analytic mode solver [13 to numerically cal- 
culate the fields transmitted through apertures in ideal 
metal sheets with cartesian symmetry. Figure [l] illus- 
trates the shadow-side longitudinal field at a distance 
A/ 10 after the metal sheet of thickness D = A/2, pro- 
duced by a left-handed or cr+ circularly-polarized plane 
wave with zero orbital angular momentum mi = inci- 
dent on a square aperture of length = Ly = X. The 
amplitude [Fig. [ija)] demonstrates electric field enhance- 
ment at the aperture surfaces. While the incident longi- 
tudinal field component has zero phase helicity = 0, 
the shadow-side phase contains a topological charge of 
rriz = +1 (clockwise, red- white-blue) [Fig. [l](b)]. The 
orthogonal (T_ or right-handed circularly-polarized field 
produces a vortex phase front with opposite topological 
charge = — 1 (clockwise, blue-white-red) [not shown]. 
Sharp aperture edges are responsible for the strong longi- 
tudinal field enhancement and large transverse field gra- 
dient. The polarization-dependent vortex phase repre- 
sents the complex exponential shown in Eq. [3j 
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FIG. 1: Numerically calculated (a) amplitude and (b) phase of 
the scattered longitudinal field produced by a left-handed or 
(7+ circularly-polarized plane wave incident on a square metal 
aperture of width Lx — Ly — X in a metal sheet with thick- 
ness A/2, (c) Poynting vector (blue) and angular momentum 
vector (red) of the scattered longitudinal field. The arrows re- 
verse direction depending on the incident orthogonal-circular 
polarization illuminating the metal sample with square aper- 
ture. 





FIG. 2: The reconstructed amplitude and phase of scattered 
longitudinal fields produced by (a-b) square and (c-d) round 
apertures in aluminum film when illuminated with circularly- 
polarized light. The incident light wavelength is 600GHz (A = 
500/im) and the aperture width and radius are approximately 
200/im. The aluminum has a thickness of 80/im. 



If we decompose the total angular momentum into con- 
tributions in the longitudinal and azimuthal directions 
< J >= r X (E X B*) + c.c. = Jzk + 7^0, then an- 
gular momentum in the z-direction is considered a 
paraxial term and the second term involves spin-orbit 
interactions described here due to the longitudinal field 
components [19 . The azimuthal component describes 
nonparaxial effects and changes sign with the incident 
circular polarization handedness. In Fig. [TJc) we illus- 
trate the Poynting vector (blue arrows) and angular mo- 
mentum vector (red arrows) of the scattered longitudinal 
field. 

We experimentally measure the scattered longitudinal 
electric fields from subwavelength apertures using a THz 
near-field electro-optic detection method and focused 
probe beam, which provides full vector characterization 
of the transmitted shadow-side electric field with 10/im 
resolution [20 . Scattering by circularly-polarized fields 
is accurately reconstructed using linearly-polarized THz 
signals. Figure [2] shows transmitted longitudinal field 
amplitudes and phases produced by incident circularly- 
polarized light at 600 GHz (A = 500/im) on circular (ra- 
dius a = 100/im) and square (length = Ly = 200/im) 
apertures. Both square [Figs. [2|a,b)] and round [Figs. 
|2|c,d)] amplitudes show field enhancement due to inter- 
action with aperture surfaces at the metal aperture edges, 
and the formation of an on-axis phase singularity. Pairs 
of opposite-sign singularities appear off-axis in Fig. |2|d), 
which we attribute to the non-normal angle of incidence 
between the signal and the sample. Moreover, these 
off-axis phase singularities associated with misalignment 
change in location depending on the incident <j+ or cr_ 
orthogonal circular polarization. 

In the longitudinal component of electromagnetic 
fields, orbital angular momentum phase singularities 
combine with the polarization-dependent phase singu- 
larities associated with spin angular momentum and in 
the scattering by subwavelength apertures, the effects of 
these phase singularities is significant. We consider input 
Laguerre- Gaussian profiles with index p = 
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where C is a normalization constant such that 
JlLG^^^^'^dA = 1, or C = (7^mz!)-^/^ and the radial 
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FIG. 3: Amplitudes and phases of the longitudinal z- 
component of transmitted fields produced by Laguerre- 
Gaussian beams of topological charge mo = +1 and (a-b) 
left-handed circular polarization and (c-d) right-handed cir- 
cular polarization. 



coordinate pra is normalized such that the mode field ra- 
dius = 1. 

Figure [3] shows the amplitude and phase of the 
scattered-field z-components after transmission through 
a square aperture given incident orthogonal circularly- 
polarized Laguerre- Gaussian beams with topological 
charge of mi = +1. The combination of left-handed 
circular polarization and the incident vortex phase pro- 
duces a net transmitted topological charge of = +2 
[Fig. |3jb)]. In contrast, the combination of right-handed 
circular polarization with similar phase vorticity cancel 
and produce a flat-phase = field [Fig. [sj^d)]. A 
comparison of the amplitudes indicate that right-handed 
circular polarization [Fig. [sj^c)] produces a shadow-side 
on-axis constructive maximum, whereas the amplitude 
of the left-handed circularly-polarized scattered field is 
strictly zero on-axis due to the on-axis phase singulari- 
ties [Fig. |3ja)]. Therefore, addition of a helical phase- 
front breaks chiral symmetry and the transmitted spatial 
beam profiles associated with orthogonal circular polar- 
izations are no longer mirror images. 

The coupling between spin polarization and orbital an- 
gular momentum phase fronts is observed in the total 
energy that is transmitted through apertures. We use 
the result of [21] to evaluate the transmission immedi- 
ately after a circular sub-wavelength aperture within an 
infinitely-thin film of metal 
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where is the angle that measures between the electric 
field vector and the x — z plane, a is the radius of the aper- 
ture, p and (f) are the cylindrical coordinates of the scat- 
tered field, and the incident angle measures between 
the z-axis and the direction of incidence, k. Equation |5] 
provides an approximation of the scattered longitudinal 
fields for the condition 27ra/A < 1. 

We decompose incident Laguerre Gaussian beams [Eq. 
|4] into plane waves with different wavevector k to solve 
for the scattering given by Eq. [5] and a± circular- 
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where the coefficient k, = ^/l — {kp/kmY contains the 
normalized wavenumber km = 2n/ X^/\ml\~^^^ and the 
cutoff transverse wavenumber is kc. For this investiga- 
tion, we consider that plasmons are excited on the inci- 
dent metal surface so that kc = oc; a cutoff wavenumber 
of kc = km implies that, for example, a dielectric coating 
prevents surface waves or plasmons from propagating on 
the incident metal surface. This k-space relation, Eq. [6j 
indicates that the scattered field amplitudes couple into 
Bessel functions of order mi and ± 2, while the longi- 
tudinal phase vorticity remains described by the relation 
mz = mi-\-ms = mi±l. The electromagnetic spin-orbit 
interaction exists in the Bessel term Jmi±2-, which indi- 
cates that the energy transmitted through the aperture 
depends on both spin angular momentum and orbital an- 
gular momentum. 

In Fig. [4] we show the difference in transmission as- 
sociated with each orthogonal polarization AT = (T+ — 
T_)/(T+ + T_) for varying orbital angular momentum 
topological charge mi as a function of aperture diam- 
eter or width = Ly = 2a, where the transmission is 

T± = J \Ei'^''''^^\'^dA. The aperture width is normalized 
by the input beam mode field diameter and since there 
is no cutoff transverse wavenumber /cc = 00, our analyti- 
cal calculation is independent of wavelength. Theoretical 
predictions are shown in Fig. [4]^a). In Fig. [4jb), we plot 
the difference in transmission from numerical simulations 
for metal sheet thickness D = X and polynomial curve- 
fit lines. The difference in transmission AT represents 

Difference in Transmission due to Polarization 
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FIG. 4: Difference in transmission between orthogonal- 
circular ly polarized vortices AT = (T+ - T_)/(T+ + T_) as a 
function of normalized hole aperture size for different incident 
topological charge mi . (a) Theoretical prediction model given 
round aperture in infinitely-thing metal sheet, (b) Numerical 
calculation for a square aperture given finite sheet thickness 
D = X. 
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the spin-orbit interaction via electromagnetic scattering 
and is therefore zero when there is no orbital angular 
momentum or = 0. 

The difference in transmission AT for a single incident 
vortex = 1 is less than 10%, while that for higher- 
order vortices exceeds 50% depending on aperture size. 
Both numerical calculations and theoretical analysis pre- 
dict similar aperture sizes for maximum spin-orbit in- 
teraction or maximum |AT|. We observe a maximum 
spin-orbit interaction for mi = 2 when aperture sizes ap- 
proach zero, and the maximum spin-orbit interactions 
occur for increasing rui at increasing aperture widths. 
Our theoretical prediction strongly underestimates the 
difference in transmission, particularly at small aperture 
sizes, and this is not reconciled by changing the metal 
sheet thickness in numerical simulations. Discrepancies 
arise from the comparison between cylindrical (theoret- 
ical) and cartesian (numerical) symmetry, however this 
point does not entirely resolve the differences described. 

Although the maximum change in the longitudinal 
field topological charge due to polarization and scat- 
tering at a single aperture is = ±1, we demonstrate 
our ability to manipulate and control spin-orbit inter- 
actions by producing shadow-side fields with arbitrary 
phase vorticity. The initial conditions that produce fields 
of arbitrary phase fronts are not uniquely defined unless 
there exists constraints placed on the amplitude of the 
scattered fields. We use a "necklace" arrangement of n 
apertures, where n corresponds to the desired phase vor- 
ticity rriz^ and prepare the input field polarization as a 
coherent superposition of azimuthally gap and radially 
grp polarized fields, or 

1 
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There is zero angular momentum in a e^^-polarized 
field that has a radially-symmetric field amplitude, al- 
though e^^ carries phase singularities that are visible 
when decomposed into the circular-polarization basis, or 
Grp {(j)) = [exp(— z0)(T+ ± exp(i^)(T_]/>/2. 

AP 

With a e^'^-polarized field centered and incident on an 
azimuthal arrangement of n apertures, each aperture pro- 
duces a single topological charge whose sign is determined 
by the local ^-handedness in Eq. [7| and the total charge 
associated with the entire scattered field is = nm^. 
In Figs. (Sja) andjsjd) we show n = 2 and n = 4 equally- 
spaced square holes separated by and with dimensions 
Ax — Ay = Lx = Ly = 2/3A. We numerically calcu- 
late the transmission assuming a metal sheet thickness 
D = A/4. Figures |5|b,e) show the amplitudes and Figs. 
[5|c,f) illustrate the phases that demonstrate singulari- 
ties rriz = -\-2ms and = +4ms. Incident fields have 
Laguerre- Gaussian magnitudes |LG(^^=i)|, although we 
observe that the amplitude of the incident field spatial 
beam profile does not substantially affect the phase of 
the scattered fields. 




FIG. 5: (a) Two apertures with dimensions and separation 
Ax = Ay — Lx — Ly = 2/3A and correspondingly scattered 
(b) amplitude and (c) phase produced via numerical calcu- 
lations. (d)Four apertures with dimensions and separation 
Ax = A^ = = = 2/3A and the corresponding scat- 
tered (e) amplitude and (f) phase produced via numerical 
calculations. The metal sheet has thickness D — A/4. 



In conclusion, electromagnetic spin-orbit interactions 
that occur via scattering are significant when vortex 
phase singularities associated with orbital angular mo- 
mentum propagate at subwavelength distances to scat- 
tering edges. Here we have demonstrated that phase sin- 
gularities associated with circular polarization and vor- 
tex phase fronts combine in the longitudinal electric field 
components. With the addition of orbital angular mo- 
mentum, the energy that is scattered through round 
and square subwavelength apertures is polarization- 
dependent. Depending on the aperture size and incident 
phase vorticity, the difference in transmission due to po- 
larization can exceed 50%. Our research indicates that 
planar asymmetric or chiral metal nanostructures that 
change the polarization of scattered fields also impart 
a phase vorticity and this explains previously observed 
polarization-dependent spatial beam profiles [7[ El [9]. 
Our results illuminate new considerations for manipulat- 
ing plasmons or surface waves, and designing and aligning 
metamaterials. 
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